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On the Binary Frequency of the Lowest Mass Members of the Pleiades with 
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ABSTRACT 

We present the results of a Hubble Space Telescope Wide Field Camera 3 imaging 
survey of 11 of the lowest mass brown dwarfs in the Pleiades known (25 — 40 M Jup ). 
These objects represent the predecessors to T dwarfs in the field. Using a semi-empirical 
binary PSF-fitting technique, we are able to probe to O'/OS (0.75 pixel), better than 2x 
the WFC3/UVIS diffraction limit. We did not find any companions to our targets. 
From extensive testing of our PSF-fitting method on simulated binaries, we compute 
detection limits which rule out companions to our targets with mass ratios of q > 0.7 
and separations a > 4 AU. Thus, our survey is the first to attain the high angular 
resolution needed to resolve brown dwarf binaries in the Pleiades at separations that 
are most common in the field population. We constrain the binary frequency over this 
range of separation and mass ratio of 25 — 40 M Jup Pleiades brown dwarfs to be <11% 
for la (<26% at 2a). This binary frequency is consistent with both younger and older 
brown dwarfs in this mass range. 


1. Introduction 

Hundreds of brown dwarfs have now been identified in the solar neighborhood through wide- 


forming regions (e.g., 

1 - ? ~ 

Epchtein et al. 

1997 

^ ^ j 

Delfosse et al. 

1997 

Chiu et al. 2006: Allers et al. 2006: 

Bihain et al. 2006: P 

„eid et al. 

2008 

Bihain et al. 

2010; 

Burningham et al. 

2010: Cushing et al. 

2011: Liu et al. 

2011 

Lodieu et al. 

2012; 

Burningham et al. 

2013). The study of brown dwarf 


binarity is a fundamental tool for testing theory, given that t he statistical pro pertie s of binaries 


2009 


probe formati o n scenarios in the very low-mass regime (e.g., iBureasser et al.1120071 : iBate 
Luhmanl 2012; Bate 2012l). For the past decade, HST and ground-based adaptive optics (AO) 
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(e.g. 

Martin et ah 

1998; Burgasser et ah 

2003: Bouv et ah 

i 

2003; 

\ ^ ~^ v ' p -- , 

Burgasser et ah 2006; Liu et ah 

2006 

and in young (1 — 10 Myr) star-forming regions such as 

Upper Sco ( 
Taurus (e.g. 

Kraus et ah 

2005; 

Bouv et ah 2006b: 

Biller et ah 2011: 

Kraus & Hillenbrand 

2012), 

Kraus et ah 

2006; 

KonoDackv et ah 2007: Todorov et al 

2010; Kraus & Hillenbrand 

2012; 

rodorov et ah 2014 

). and 

Chamaeleon I (e.g. 

Neuhauser et ah 2 

002; 

Luhman 2004: Lafreniere et ah 

2008; 

Ahmic et ah 

2007; 


Luhmanll2007l ). Multiplicity studies ha ve also been performed in older (~400 Myr) regi ons such as 
Coma Ber, Praesepe, and the Hyades (jKraus fc Hillenbrandl 120071 : iDuchene et al.l 1201. 'll ). 


Previous work has shown that t he b inary frequ ency decreases and typical mass ratios increase 
going to lower mass primaries (IBurgasser et ah 20071) • One surprising finding is that these properties 
apparently differ between young and old binaries, with the binary frequency enhanced at young 


ages by a factor of ~2x (e.g., [Lafreniere et al. 20081) and with wide separations (~10 — 1000 AU) 


being much more common as compared to field brow n dwarf binaries that are rarely wider than 
10 AU (e.g. IBurgasser et al.l 120061 : IClose et al.ll2007l ). An unambiguous physical explan ation for 
this diff erence is lack ing, as even relatively wide binaries in young star-forming regions (iLuhman 
2004j; [Lu hman et ah 20091 ) are not expected to incur dynamical interactions of sufficient intensity 
to reduce their frequency and truncate their separation distribution. 

The Pleiades open cluster serves as an important bridge between the youngest (1 — 10 Myr) 


age of »125 Myr (Stauffer et ah 

1998; Barrado v Navascues et ah 2004) and distance of 136. 

2=1= 1.2 

’ies in 

pc (Melis et al 

. 2014|b There are manv surveys that have searched for brown dwarf binai 

the Pleiades 1 

dartin et ah 2000 

Dobbie et ah 

2002; Jameson et ah 2002: Nagashima et ah 

2003; 

Moraux et ah 

2003: Bouv et ah 

2006a). However, there are only 4 Pleiades brown dwarfs with 


primary masses < 40 M Jup that have been searched for companions to date (|Moraux et al.1120031 : 
Bouv et a h I2006al ). At such masses, these objects will cool to T dwarfs at ages of the field popula¬ 
tion. 


In this work, we triple the number of low mass Pleiades brown dwarfs searched for companions, 
surveying a sample of 11 previously unobserved L dwarfs in the Pleiades using HST/WFC3. We 
computed detection limits for our sample using a binary fitting technique and Tiny Tim PSF 
models. We compared our binary frequency to the observed frequencies for brown dwarfs at similar 
masses in Taurus, Chamaeleon I, Upper Scorpius, and the field population. 


2. Observations 
2.1. Sample 

We obtained images of 11 Pleiades brown dwarfs using the Hubble Space Telescope ( HST) with 
the UVIS channel of Wide Field Camera 3 (WFC3/UVIS) in January and February of 2012 (GO 
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12563, PI Dupuy). Our sample consists of the faintest (K > 16 mag), la test type (>M 9) members 
of the Pleiades known in earl y 2011. According t o BT-Settl models of lAllardl (120141 ) tied to the 
COND evolutionary models of Parade et ah (1997, 1 998! . 2Q03j), the estimated masses of our sample 
are 25 — 40 M Jup based on their it'—band magnitudes and the age of the Pleiades. When defining 
our sample, we considered objects bona-fide members of the Pleiades if they had proper motion 
indicating cluster membership and spectra with low surface gravity features or lithium absorption. 
Our sample is listed in Table [H along wi th 4 targets from pre vio us HST/ACS and H ST/ WFPC2 
observations of Pleiades brown dwarfs by iMartfn et al.1 (|20031 ) and iBouv et al.1 (|2006af ) that match 


(Bihain et al. 

2006; 

Casewell et al. 

2007: 

Lodieu et al. 

2012) 

and BRB 29 have spectral types L0-L4.5 from 

Bihain et al. 


(i20ld ). 


2.2. HST/ WFC3 Imaging 


We obtained 2 exposures each in filters F814W and F850LP for each target star. One image 
of brown dwarf BRB 17 was lost due to a pointing error so we had a total of 43 images. The target 
stars are positioned near the center of the full field of view at ~250 pixels from the bottom of chip 
1. We chose a longer exposure time of 900 s in F814W filter, where we are sensitive to tighter brown 
dwarf binaries because of the smaller PSF. We also obtained 340 s exposures in F850LP to confirm 
the presence of any candidate companions and measure their colors. The full width half maximum 
of the PSF is «1.84 pixels in F814W and rT. 96 pixels in F850LP according to the WFC3 data 
handbook 

We inspected each image for cosmic rays hits, identified as rays or streaks with high counts 
but not resembling WFC3 point sources. We found 6 of the 43 images had cosmic ray hits within 
5 pixels of th e target star. We use the Laplacian Cosmic Ray Identification algorithm LAC0SMIC 
(jvan D okku m 200 ll) to remove cosmic rays from a 200x200 pixel area on the detector centered 
on the target star. LAC0SMIC replaces each pixel with the median of the surrounding pixels in 
an iterative procedure. Visual inspection after the fact confirms that we successfully cleared all 
obvious cosmic ray hits except for a single image of brown dwarf BRB 23 in F850LP due to a 
cosmic ray hit through the center of the peak of the target. We excluded this image of BRB 23 in 
the subsequent data analysis, therefore leaving us with 42 images total for the rest of our analysis. 

We computed aperture photometry of our targets from the pipeline calibrated, geometrically- 
corrected, dither-combined (drz) images. We calculated our aperture photometry using the APER 
task from the IDL Astronomy User’s LibrarjH for an aperture radius of 074 and a sky annulus of 
074 — 078. We converted the flux in our aperture to a Vega magnitude using zeropoints of 24.57 mag 


: http: //www. stsci . edu/hst/wf c3/documents/handbooks/currentIHB/ 
■'http: //idlastro. gsf c .nasa.gov/homepage . html 
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for the F814W filter and 23.20 mag for the F850LP filter provided in the HST/ WFC3 webpagejf). 
To determine our photometric uncertainties, we first constructed an error image for each image, 
accounting for read noise and poisson noise. Using a Monte Carlo approach, we determined our 
photometric errors from 10 4 iterations of the APER task after adding random Gaussian noise to the 
image in each iteration. The resulting F814W and F850LP photometry for our targets is listed in 
Table [2j 


3. Image Analysis 

3.1. Point Spread Function (PSF) Model of WFC3/UVIS 


In order to search for close companions to our targets, we began by fitting a model Tiny 
Tim (Krist et aJj 2011) point spread function to our imaging data. To create the most accurate 
model we specified the exact coordin ates of our target and used an input spectrum of 2MASS 
J00361617+1821104 ( Reid et al. 2000 . L3.5). We set the defocus parameter in Tiny Tim to the 
model defocus provided on the Space Telescope Science Institute webpageQ for each image of each 
target. The model defocus is computed to account for breathing, according to the telescope tem¬ 
perature data. 

We sampled the Tiny Tim PSF at 5x the pixel scale ((y/04 pixel -1 ) of WFC3/UVIS1. To 
simulate sub-pixel shifts of our targets we bilinearly interpolated to an arbitrary fractional pixel and 
then binned down to pixel scale of WFC3. We used the Nelder-Mead downhill simplex method from 
Press! (19881). which is the AMOEBA algorithm in IDL, to minimize the y 2 , varying the (x, y ) position 
and flux normalization until finding the best fit. We computed y 2 as ((image-model)/noise) 2 , where 
“noise” is the noise image provided by the WFC3 reduction pipeline. We ran the AMOEBA algorithm 
tw ice, sta rting the second run at the end point in parameter space of the first run, as recommended 
by Press (1988). We fit a ±10 pixel cutout region centered on the target star. 


We found average residuals after subtracting the best-fit Tiny Tim model of 5% and 6% 
for F814W and F850LP images, respectively. We computed residuals of our fits as the average 
fractional offset between the image and the model. The majority of the residual flux using the 
Tiny Tim model was at instrumental position angles of 30 — 50° and 150 — 170° in both the F814W 
and F850LP filters (Figure [T]). If we searched for faint companions using the TinyTim PSF model 
and our binary fitting technique detailed below, we found that this systematic residual flux led to 
spurious detections of companions at these position angles. 


Therefore, we instead computed a single optimal semi-empirical PSF model that minimized 
the residuals across all images by modifying the Tiny Tim model. We iteratively solved for a 5x 


a http://www.stsci.edu/hst/wfc3/phot_zp_lbn 
a http://www.stsci.edu/hst/observatory/focus 
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over-sampled additional component image to be added to the Tiny Tim model. The best guess of 
this additional component at each pixel was computed as the median across all normalized images 
of the data minus the previous iteration’s PSF model. We computed a semi-empirical PSF model 
as the Tiny Tim model at the mean position of our targets with this additional component added 
in. 


Using our semi-empirical PSF model, the final residuals of our fits were improved by 5x to 
«0.9% and «2.3% for F814W and F850LP, respectively (Figured]). Most importantly, we no longer 
see the concentrated residual flux at position angles of 150 — 170°. We use our semi-empirical PSF 
model in all subsequent analysis. 

The method of fitting binaries is the same as described above, but instead of using a single 
model we use two co-added models. As before, the AMOEBA algorithm minimizes the y 2 between 
the image and co-added semi-empirical PSFs. We varied six binary parameters: the primary’s 
position on detector, the flux normalization between the primary star and the PSF model primary, 
the binary separation, the position angle, and flux ratio between the primary and secondary. 


3.2. Quantifying False Positives 

If we run our binary fitting code on a image of a single star, we recover binary parameters of 
false positive companions. By definition, these detections reveal the distribution in separation and 
flux ratio of the false positives we would find while searching for companions in our imaging data. 
To characterize the false positives for our WFC3 data, we fit images of our target stars using our 
binary fitting technique from H3.il We scale all images to either the median or minimum S/N of 
our sample by adding in gaussian noise (Table [3]). This allows us to put our sample on a common 
scale for our simulations. For each target star, we start with 150 random initial guesses, uniformly 
distributed in (x,y) from 0.1 — 5 pixels, and flux ratios from 0 — 5 mag. 

We show the resulting distribution of separations vs flux ratios of recovered false positives in 
Figure [2j The brightness of false positives increases with decreasing separation. At the tightest 
separations (<0 , /02, <0.6 pixel), we find that near unity flux ratio false positives are the most 
common. At wider separations (>1.5 pixels, >(X , 06), we find that almost all false positives are 
found with large flux ratios of 3 — 5 mag. This is expected, as the binary fitting code is required to 
return a position and flux normalization for a secondary even if one doesn’t exist. In other words, 
the single WFC3 PSF can be fit with a model of a high flux PSF and a very low flux PSF added 
in to fit any small leftover residuals. 
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3.3. Artificial Binary Simulations 

In order to compute detection limits for our survey, we generated artificial binaries at random 
separations of 0.3 — 5 pixels (0 / . / 018 — O'/2), position angles of 0 — 360°, and flux ratios of 0 — 5 
mag. We created these artificial binaries by shifting, scaling and co-adding randomly selected pairs 
of actual images together. Given that the marginally sampled WFC3 PSF (FWHM <2 pixels) 
hinders the accuracy of linear interpolation at sub-pixel shifts, we shift the secondary star relative 
to the primary star in integer pixel steps. We scaled the image of every primary to a common 
S/N by adding noise, thus degrading the image to lower S/N. We scaled the secondary to a S/N 
appropriate for the randomly chosen flux ratio of the artificial binary. 

Given the integer pixel shifts, there are fixed separations and position angles allowed by the 
possible image pairings. These integer pixel shifts can result in non-integer artificial binary separa¬ 
tions because the sub-pixel position for each image varies. Out of all possible pairings we selected 
a subset of 4800 artificial binaries that are distributed uniformly in log separation, flux ratio, and 
position angle. We ran two sets of simulations for each filter, scaling primaries alternatively to the 
median S/N and the minimum S/N of our images (Tabled]). Only half the images were used for 
the median S/N simulations, given that we only scaled images down in S/N, never up. 

We then blindly fitted for the binary parameters of our artificial binaries using a double PSF 
model as described in T3.ll using 150 random initial guesses. The best-fit values for each parameter 
are calculated as the mean of the resulting 150 runs of our binary fitting code parameters where 
runs with outlier y 2 were excluded from the average. 


3.4. Deriving False Positive Curves 

The binary parameters recovered in our artificial binary simulations contain a mix of both 
detections and false positives. To assess the likelihood of a given binary fit being a detection, we 
compared our distribution of false positives from 113.21 and our fits to artificial binaries from H3.3I 
to measure our false positive curve, i.e the largest flux ratio before the recovered secondary star 
becomes indistinguishable from a false positive at a given separation. 

We considered the artificial binaries and false positives in a given separation and flux ratio 
range, using 0.1 dex pixel bin widths and 0.3 mag flux ratio bin widths, respectively. In each 
separation bin we normalized the histogram of false positive flux ratios to the histogram of recovered 
artificial binary flux ratios by conservatively assuming that any artificial binaries with recovered 
flux ratios larger than the median false positive flux ratio Am cr ; t were most likely false positives 
themselves. We computed this normalization factor as Q 5 ” , where rif p is the total number of false 

positives and n is the number of artificial binaries with flux ratios >Am cr it. After normalization, we 
computed the false positive fraction as a function of flux ratio as 1 — We repeat the procedure 
above for each separation bin. This procedure is depicted in Figure [3] for the 0.79 — 1.0 pixel 
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separation bin. 

With the procedure detailed above, we computed false positive curves at the median and 
minimum S/N of our images for the F814W and F850LP filter as shown in Figure 01 Each of our 
false positive curves are representative of a single , S/N given that we scale our all our images to a 
common S/N for each set of simulations. 


3.5. Deriving Contrast Curves 


We computed contrast curves that correspond to the largest flux ratio companion that our 
binary PSF fitting technique can recover accurately at a given separation. A binary is considered 
“recovered” if the best fit parameters are within 0.2 pixels and 1 mag of the input (x, y ) positions 
and flux ratio, respectively. We binned our simulated binaries by separation and flux ratio with 
bin widths of 0.1 dex pixels and 0.3 mag, respectively. In each bin, we computed the completeness 
fraction as the number of artificial binaries that are recovered divided by the total number of 
artificial binaries in the bin. We define our contrast curves as the flux ratio bin at a given separation 
where the completeness fraction is 90% determined by the interpolation of the binned results. We 
computed contrast curves at the median and minimum S/N of our targets (Table [3]) for the both 
F814W and F850LP filters. 


Figure [5] shows our resulting contrast curves. We are able to recover tight (<0 / /04, <1 pixel) 
binaries with flux ratios <1 mag. At wider separations we recover binaries 3 — 5 magnitudes fainter. 
We also constructed a contrast curve with a stricter recovery requirement to be within 0.3 mag of 
the input. This leads to a contrast curve that reaches in to binary separations of 0 , . , 035 (0.9 pixels) 
and is identical to our default recovery requirements outside Of'055 (1.4 pixels). A flux ratio of <1 
mag for our targets corresponds to a mass ratio q > 7 which allows us to rule out the possibility 


of Pleiades brown dw arf binaries si mil ar to held brown dwarf binaries, since the latter mostly have 
q ~ 1 (see review by IBur gasser et al. 2007]). This means that a stricter flux ratio requirement of 
<0.3 mag for constructing our contrast curves is unnecessary. Thus, our PSF fitting technique is 
able to recover artificial binaries as tight as 0 , . / 03, well inside the diffraction limit (« ^ X/D ) 


Given that each target in our sample has a different S/N, we interpolated over the measured 
median and minimum S/N curves to compute a contrast curve for each target. We conservatively 
fixed the contrast curve for our targets with S/N higher than the median S/N to the median 
S/N contrast curve. Our detection limits in F814W and F850LP mag for each target are shown 
in Table 0J These detection limits have lower contrast and are more conservative than the false 
positive curves, as expected. Finally, we convert our cont rast curves from F814W and F850LP 
m agnitudes to masses using BT-Settl models I Allard] (l2014l| tied to the COND evolution m odels 
of iBaraffe et al.l ( 2003h . We assumed an age of 125 My r (jBarrado v Navascues et al.l 1 20041 ) and 
distance to the Pleiades of 136.2 pc (jMelis et al. [2014]). Figure 6 shows the 90% completeness 
contrast curve for each target as a function of mass ratio ( q ) and projected separation (a) in AU. 























We use only the F814W contrast curve for our constraint on the binary frequency due to higher 
S/N, larger contrast, and closer limiting separation than our F850LP contrast curve. 


3.6. Completeness Maps 

Similar to how we derive contrast curves in 33.51 we derive a median and a minimum S/N 
completeness map for the F814W and F850LP filters. Each completeness map represents the 
probability that a companion with a given separation and flux ratio would have been detected 
(Figure [7]). The procedure for deriving completeness maps is exactly as deriving a contrast curve in 
113.51 except that we compute the completeness fraction at every separation and flux ratio bin. We 
computed a completeness map for each target similar to 33.51 by interpolating over the median and 
minimum S/N completeness maps. We conservatively fixed the completeness maps for our targets 
with S/N higher than the median S/N to the median S/N completeness map. Our completeness 
maps for several targets are shown in Figure [TJ 


4. Results 


4.1. L Dwarf Binary Frequency of the Pleiades 


We found no companions in surveying 11 brown dwarf members of the Pleiades with K > 16.0 
mag. Our F814W contrast curves demonstrate that we could have detected companions with mass 
ratios of q > 0.5 at separations a > 10 AU and q > 0.8 at a > 4 AU (Fi gure El). Most k nown very 
low mass binaries are sharply peaked towards mass ratios q ~ 1 (jBurgasser et al.l 120061 ; ILiu et al 


2010 ). Furthermore, our detection li mits probe do wn to separations a ~ 4 AU, near the peak of 


the observed binary distribution (jBurgasser et al.ll2006l h Thus, our detection limit s are sensitive to 


the majority of binaries expected from the ob s erved field populatio n of T dwarfs (jBurgasser et al 


20031 . 120061 : lOelino et al.l 1201 ll : ILiu et al.ll2012l : iRadigan et al.ll2013l ). 


We estimated the binary frequency for the Pleiades by comparing our completeness maps ( 33.6p 
to various random simulated populations of binaries. Each population of binaries had an adopted 
eccentricity, mass ratio and separation distribution, with semi-major axes of <25 AU in accordance 



longitude of ascending node, mean anomaly, and argument of periapsis, and an a sin?) distribution 
for inclination. We projected each binary on sky from the population with 10° randomly chosen 
orbits. We compared each of these 10 5 orbits to each completeness map of each target. The 
probability for detecting a binary was given by our completeness fraction at the separation and 
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mass ratio of the binary from the completeness maps (Figure [7|. We averaged over all probabilities 
and computed a single a verage probabili ty Udetect ability”) to recover a companion for each target 
star (Table[6]). Similar to Aberasturi et ah (120141 1. we then summed over these average probabilities, 
and found that if all our targets had companions we should have detected 7.6 binaries for the log 
normal distribution o f se mi- major axes . We also used a linear (flat) semi-major axis distribution 
to be consistent with Aberasturi et ah (120141 ). finding virtually no difference in the total number 
of binaries we should have detected (8.1). The lack of detections implies a binary frequency upper 
limit of <11% for l a (<26% at 2a) using the recommended Jeffrey’s distribution for small n 


(jBrown et al. Il200lh . lAberasturi et al.1 (|2014T ) computed a binary frequency for >T5 primaries in 


the solar neighborhood of <16%-<25% using the Clopper-Pearson interval at 95% confidence using 
the same log normal and uniform separation distributions. This is comparable to our own binary 
frequency upper limit of <26% at 2a (« 95% confidence). 


4.2. Binary Frequency vs Age for Wide (>10 AU) Companions 


According to the evolution models of lBaraffe et, al.1 (20031), our sample of Pleiades L dwarfs are 


expected to evolve to T e g = 700 — 1300 K (i.e., T0-T8 spectral types) at ages of 0.5 — 5.0 Gyr. At 
younger ages of 1 — 10 Myr, our sample would have had temperatures of 2300—2750 K (i.e. M7—M9). 


in Taurus (Todorov et ah 2014: 

Kraus & Hillenbrand 2012; Kraus et ah 

2006; 

KonoDackv et ah 

2007; Todorov et ah 

20ld), Chamaeleon I ( 

Luhrnan 

2004; 

Lafreniere et ah 

2008; 

4hmic et ah 

2007; 

Luhrnan 

2007; 

Neuhauser et ah 

2002), Upper Sco 

(Biller et ah 2011: Kraus & Hillenbrand 

2012) 


and the field ( Burgasser et ah 20061 ). Taurus. Chamaeleon I and Upper Sco are regions with objects 


all at the same distance, thus aiding the comparison. 


It is possible that the different cluster stellar densities in which brown dwarfs form could 
affect the binary frequency, hinderin g a direct compari son between field and young brown binary 
frequencies as done here. However, King et ah (120121 ) find that the binary frequency for stars 
with masses of 0.1 — 0.3 Mq did not vary measurably over nearly 20 x in density for five young 
regions (Taurus, Chamaeleon I, Ophiucus, IC 348, and the Orion Nebula Cluster). Figure [8] and 
Table [5] summarizes these comparisons of the binary frequency at different ages. In contrast to o ur 
estimate of the binary frequency in H4.ll here we used only the methods of Burgasser et ah ( 20031 ) 
for computing the binary frequency of these different clusters and the held in order to keep the 
statistical analysis the same. 


For constraining our binary frequency of Pleiades at wider separa tions a > 10 AIK 4 brow n 


dwarfs observed by the I75T/WFPC2 and HST/ACS surveys of Martin et al.1 (|2003l ) and lBouv et al 


(|2006al ) were combined with our own observations for a larger sample size of 15 objects. These 
4 brown dwarfs match our K > 16.0 mag cutoff and conservative Pleiades cluster membership 


criteria, i.e. that the target must have proper motion indicating cluster membership and a spectral 
type >M9 (see H2.1|h Brown dwarfs PLIZ 28 and PLIZ 2141 were observed with HST/ACS by 
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Bouv et al. (2006a) with detection limits that ruled out companions for mass ratios q > 0.45 at 
sep arations a > 7 —12 AU. Brown dwarfs Roque 30 and Roque 33 were observed with HST/ WFPC2 
by [Martin et al. (|2003l ) and similarly they ruled out companions for mass ratios q > 0.5 and 
separations a > 10 AU. The HST /ACS and HST/ WFPC2 observations have comparable detection 
limits to our own detection limits of q > 0.6 at separations a > 10 AU. Thus, with a combined 
sample size of 15 low mass Pleiades brown dwarfs and no binaries detected, we computed an upper 
limit on the binary frequency of <7.0% (lu) for mass ratios q > 0.6 and separations a > 10 AU. 

The sa mple of you ng b rown dwarfs observed by HST/ WFPC2 and AO surveys (see Table [5]) 
compiled in iTodorov et al. (1 20141 ) and references therein includes all targets with spectral types 
>M4. The detection limits for these surveys are generally sensitive to companions with separations 
a > 10 AU. In an attempt to constra in the masses of the primaries to <40 M Jup , we included only 
primaries in the Todorov et al. (2014) sample with spectral types >M7 (see Table [5]). Note that for 
young (<10 Myr) brown dwarfs mass estimates at young ages are still uncertain and could have large 
uncertainties due the lack of a well measured T e g scale for these stars and uncertain atmospheric and 
stellar evolution models. This spectral type cut off corresponds to a mass estimate of <40 M Jup at 
age s ~ 1 Myr and m2 — 3 Myr for the Taurus and Chamaeleon I regions, respectively, according to 


the lBaraffe et al. (120031 ) models. Over this range there are 3 out of 37 binaries in Taurus and 1 out 
of 22 binaries in Chamaeleon I, which corresponds to binary frequencies of 0.0 —6.0% and 0.0—10.0% 
(lcr) respectively. We find our binary frequency upper limit of <7.0% is in agreement with binary 
frequencies for both Taurus and Chamaeleon I. One caveat is we included ca ndidate compan i ons in 
Taurus 2MASS J04414489+2301513 and 2MASS J04221332+1934392 from ITodorov et al] ( 20141 ) 
in the binary frequency computed here. If those objects are not binaries, the binary frequency of 
Taurus would be even lower (0.0 — 6.0%), still in agreement with our own binary frequency limit. 


Kraus fc Hillenbrandl (120121 ) and iBiller et al.l (120111 1 observed 10 and 18 members of Upper 


Sco with spectral types >M7 respectively and were sensitive t o companions with separations >10 
A U. Given an age of 1 1 Myr for Upper Sco dPecaut et al.ll2012ll and the spectral type-T e g relation 
of P ecau t &; Mama jeld (120131 ). >M7 spectral types correspond to <2650 K and thereby masses of 
<40 M Jup . This is comparable to our own mass range of 25 — 40 M Jup . Both previous surveys have 
detection limits q > 0.8 at separations a > 10 AU with no binaries detected. Using this combined 
sample, we estimated a binary frequency of 0.0 — 4.0% for Upper Sco, which is consistent to our 
own binary frequency upper limit of <7.0% for the Pleiades. 


Burgasser et al.l ( 20061 ) resolved 5 T dwarf binaries with separations of a = 1.8 — 5.0 AU 


out of 22 stars observed with HST/ NICMOS. They computed a Malmquist bias-corrected binary 
frequency of 8 — 19% for mass ratios q > 0.6 and separations a > 2 AU. However, to directly 
compare to our detection limits, we recomputed their Malmquist bias-corrected binary fraction 
and considered only the 2 T dwarf binaries which have projected separations of >10 AU, which 
gives a binary frequency of <3.0% for 0 binaries detected out of 17 objects observed. 


Batel (2 0121 ) performed hydrodynamic simulations of star formation that produced 27 objects 
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with masses <70 M Jup , with none ending up as binaries. Bate! (2012) quoted a binary frequency 
of 0.0 ± 5% for the mass range of 30 — 70 M Jup and a binary frequency of <7% for the mass 
range 10 — 30 M Jup . These predictions are in good agreement with our observed binary frequency 
constraint of <7.0% for separations >10 AU. 


5. Summary 

The measurement of the brown dwarf binary frequency at different ages is fundamental tool for 
testing theory, given that the statistical properties of binaries probe formation scenarios in the very 
low-mass regime. In this work, we tripled the number low-mass Pleiades brown dwarfs searched for 
companions, surveying a sample of 11 previously unobserved L dwarfs in the Pleiades, predecessors 
to T dwarfs in the field, using HST/ WFC3. We have constrained the binary frequency in Pleiades 
for the lowest known mass (25 — 40 M Jup ) and latest known type (>M9) brown dwarfs to <11% at 
la (<26% at 2cr) confidence for companions as close as « 4 AU, finding no binaries. Our survey is 
the first to probe down to separations of 4 AU at such young ages. 

Furthermore, we find our binary frequency constraints are in good agreement with observed 
binary frequencies of young star forming regions Taurus (0.0 — 6.0%), Chamaeleon I (0.0 — 10.0%), 
and Upper Sco (0.0 — 4.0%) for objects with similar primary masses of <40 M Jup , at la with 
projected separations >10 AU. Overall, our observations of the Pleiades support the evidence that 
T dwarf binaries are likely uncommon, and consistent with having the same frequency at both 
young (1 — 10 Myr), intermediate («120 Myr) and old (>1 Gyr) ages. 
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Table 1. Pleiades Sample 


Name® 

R. A. 

Decl. 

Mass 6 

K 

SpT SpT 

P.M. 


J2000.0 

J2000.0 

Mj up 

mag 

Ref 

Ref 


BRB 17 

03 

54 

07.98 

+23 

54 

27.9 

43 

16.03 

± 

0.03 

L0 

1 

2 

NPNPL2 

03 

46 

34.26 

+23 

50 

03.7 

41 

16.09 

± 

0.03 



3 

PLIZ 31 

03 

51 

47.65 

+24 

39 

59.2 

40 

16.09 

± 

0.03 



3,4 

BRB 21 

03 

54 

10.27 

+23 

41 

40.2 

31 

16.39 

± 

0.04 

L3 

1 

2 

PLIZ 35 

03 

52 

39.16 

+24 

46 

29.5 

31 

16.51 

± 

0.04 

L2 

1 

2 

BRB 23 

03 

50 

39.54 

+25 

02 

54.7 

30 

16.56 

± 

0.04 

L3.5 

1 

2 

PLIZ 161 

03 

51 

29.47 

+24 

00 

37.3 

28 

16.70 

± 

0.05 



3 

UGCS J0348+2550 6 

03 

48 

15.63 

+25 

50 

08.9 

28 

16.73 

± 

0.05 

L3+1 

8 

3,7 

BRB 28 

03 

52 

54.90 

+24 

37 

18.2 

26 

16.92 

± 

0.06 



2 

PLIZ 1262 

03 

44 

27.27 

+25 

44 

42.0 

26 

16.95 

± 

0.07 



2,4 

BRB 29 

03 

54 

01.43 

+23 

49 

57.7 

25 

17.00 

± 

0.07 

L4.5 

1 

2 

Roque 33 c 

03 

48 

49.03 

+24 

20 

25.4 

41 

16.06 

± 

0.03 

M9.5 

6 

5 

Roque 30 c 

03 

50 

16.09 

+24 

08 

34.7 

40 

16.08 

± 

0.03 



3 

PLIZ 28 d 

03 

54 

14.03 

+23 

17 

51.4 

35 

16.14 

± 

0.03 

LO.O 

1 

2 

PLIZ 2141 d 

03 

44 

31.29 

+25 

35 

14.4 

28 

16.69 

± 

0.04 



2 


Note. — “To search thes e targets by name in Simbad, add the string “Cl* Melotte 22” 


b Masses are estimated from iBaraffe et al. ( 2003) 
c Observed with HST/ WF PC2 Martin et al. ( 2003 ) 


Observed with HST/AGS [Bouv et ali (12006al) 


e UGCS J034815.64 +255008.9 


References. (1) Bihain et al. (2010); (2) 

Bihain et al. (2006 

; (3) 

Casewell et al. (20071); (5) 

Stauffer et al. (2007 

M6) 

Martin et al. 

20001 

2014a): 18) Zapatero Osorio et al. ( 

2014b) 



Lodien et all (2012 ); (4) 


f7) IZapatero Osorio et al 
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Table 2. HST/WFC3 Photometry 


Our Targets 

F814W 

(mag) 

F850LP 

(mag) 

BRB 17 

20.419 ±0.007 

19.415 ±0.011 

NPNPL2 

20.685 ± 0.006 

19.448 ±0.011 

PLIZ 31 

20.701 ±0.006 

19.524 ±0.013 

BRB 21 

21.344 ±0.010 

20.204 ± 0.023 

PLIZ 35 

21.315 ±0.010 

20.096 ±0.021 

BRB 23 

21.604 ±0.012 

20.431 ± 0.029 

PLIZ 161 

21.804 ±0.014 

20.678 ± 0.034 

UGCS J0348+2550 

21.866 ±0.015 

20.706 ± 0.035 

BRB 28 

22.177 ±0.019 

20.860 ± 0.040 

PLIZ 1262 

22.211 ±0.020 

21.086 ± 0.049 

BRB 29 

22.231 ±0.021 

21.042 ±0.048 


Table 3. Binary Simulations 


Simulation 

Filter 

S/N 

Number of 

Artificial Binaries 

Median S/N 

F814W 

93.5 

4800 

Min S/N 

F814W 

61.1 

4800 

Median S/N 

F850LP 

49.1 

4800 

Min S/N 

F850LP 

33.0 

4800 










Table 4. Detection Limits 


Target 

0"025 

0"032 

O'/040 

0"050 

0'. / 063 

0"080 

0"100 

O'.'126 

O'.'159 

O'.'180 

F814W (mag) 

BRB17 

20.66 

21.26 

21.88 

22.20 

22.48 

22.75 

22.77 

23.05 

23.69 

23.06 

NPNPL2 

20.93 

21.52 

22.15 

22.46 

22.75 

23.02 

23.04 

23.31 

23.95 

23.33 

PLIZ31 

20.94 

21.54 

22.17 

22.48 

22.76 

23.03 

23.05 

23.33 

23.97 

23.35 

BRB21 

21.59 

22.18 

22.81 

23.12 

23.41 

23.67 

23.70 

23.97 

24.61 

23.99 

PLIZ35 

21.56 

22.15 

22.78 

23.09 

23.38 

23.65 

23.67 

23.94 

24.58 

23.96 

BRB23 

21.85 

22.44 

23.07 

23.38 

23.67 

23.93 

23.96 

24.23 

24.87 

24.25 

PLIZ161 


22.64 

23.27 

23.26 

23.54 

24.13 

23.86 

24.43 

24.72 

24.45 

UGCSJ0348+2550 


22.70 

23.33 

23.32 

23.61 

23.93 

23.92 

24.49 

24.78 

24.51 

BRB28 


22.70 

23.30 

23.63 

23.92 

24.25 

24.23 

24.53 

24.76 

24.55 

PLIZ1262 


22.73 

23.34 

23.66 

23.66 

24.28 

24.27 

24.57 

24.80 

24.58 

BRB29 


22.75 

23.36 

23.68 

23.68 

24.30 

24.29 

24.59 

24.82 

25.44 




F850LP (mag) 






BRB17 



20.26 

20.55 

20.57 

20.90 

21.49 

21.14 

21.77 

21.76 

NPNPL2 



20.29 

20.59 

20.60 

20.93 

21.52 

21.17 

21.81 

21.80 

PLIZ31 



20.37 

20.66 

20.68 

21.01 

21.60 

21.25 

21.88 

21.87 

BRB21 



21.05 

21.34 

21.36 

21.69 

22.28 

21.93 

22.56 

22.55 

PLIZ35 



20.94 

21.23 

21.25 

21.58 

22.17 

21.82 

22.46 

22.44 

BRB23 



21.28 

21.57 

21.59 

21.91 

22.51 

22.16 

22.79 

22.78 

PLIZ161 



20.73 

20.92 

21.83 

22.16 

22.12 

22.40 

22.44 

22.16 

UGCS J0348+2550 



20.76 

20.95 

21.86 

22.19 

22.15 

22.43 

22.46 

22.19 

BRB28 




20.91 

21.71 

22.02 

22.31 

22.30 

22.30 

22.34 

PLIZ1262 




21.14 

21.94 

22.24 

22.53 

22.53 

22.53 

22.57 

BRB29 




21.10 

21.90 

21.88 

22.49 

22.48 

22.48 

22.53 
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Table 5. Binary Frequency vs. Age for Wide (>10) AU Companions 


Region 

Age 

Age 

Ref 

Sample 

Ref 

Aobj 

Abin 

Bin Freq 

% 

q 

Taurus 

1 Myr 

15 

1,2,3,4,5 

37 

3 

0.0 - 6.0 

>0.7 

Chameleon I 

2-3 Myr 

16 

4,5,6,7,8,9,10 

22 

1 

0.0-10.0 

>0.7 

Upper Sco 

11 Myr 

17 

2,11 

28 

0 

0.0 - 4.0 

>0.8 

This work + lit 

125 Myr 

18 

12,13 

15 

0 

<7.0 

>0.6 

Field 

0.5-5.0 Gyr 

19 

14 

17 

0 

<3.0 

>0.6 


Note. — Faint companions to brown dwarfs with separations and mass ratios greater 
than given in table are ruled out by the given detection limits for primaries with masses 
<40M Jup and separati ons > 10 AU. 


References. 


11 ) 


Kraus et al.1 ( 2006) (4) 


(2004); (7) Lafreniere et al 


( 10 ) 

(13) 


odorov et al. 


(120141): (2) Kraus fe Hillenbrandl (120120: (3) 


Konopackv et all ( 2007) (5) Todorov et. al.l ( 201C); (6) 


Luhman et al 


M 


Ahmic et al. 


Neuhauser et al. 

( 

2002 

); (ll) 

Biller et al. 

Bouvet a 

( 

2006a 

); (14) 

Burgasser et al. 


(20071): (9) Luhman (2007): 


Mill); (12) 
(2Q0®; (1® 


Martin 


et al. 


Luhman 


_ (120031 ): 

Luhman i2f)()i]): (16) 


( 2010l ): (17) IPecaut et al.1 ( 2012); (18) Barrado v Navascues et al. ( 2004 ); 


(19) A ssumed age for field T dwarfs by iBurgasser et al.1 (|2006l ) from iReid fc Hawlev 

( 120001 ). 
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Table 6. Companion Detectability 


Name 

Detectability 
log normal a 

BRB17 

70.8% 

BRB21 

68.7% 

BRB23 

69.9% 

BRB28 

67.2% 

BRB29 

66.7% 

NPNPL2 

71.4% 

PLIZ1262 

66.9% 

PLIZ161 

68.8% 

PLIZ31 

71.3% 

PLIZ35 

68.6% 

UGCS J0348+2550 

68.8% 

Total Expected Binaries 

7.6 

Binary Frequency® 

<11% 


Note. — a Binary frequency with la 
using the Jeffrey interval reco mmended for 


low n by Brown et ah (2001). 









-17- 


Bate, M. R. 2009, MNRAS, 392, 590 
Bate, M. R. 2012, MNRAS, 419, 3115 

Bihain, G., Rebolo, R., Bejar, V. J. S., et al. 2006, A&A, 458, 805 

Bihain, G., Rebolo, R., Zapatero Osorio, M. R., Bejar, V. J. S., & Caballero, J. A. 2010, A&A, 
519, A93 

Biller, B., Allers, K., Liu, M., Close, L. M., & Dupuy, T. 2011, ApJ, 730, 39 
Bouy, H., Brandner, W., Martin, E. L., et al. 2003, AJ, 126, 1526 
Bouy, H., Moraux, E., Bouvier, J., et al. 2006, ApJ, 637, 1056 
Bouy, H., Martin, E. L., Brandner, W., et al. 2006, A&A, 451, 177 

Brown, LD, Cat, TT and DasGupta, A (2001). Interval Estimation for a Binomial proportion. 
Statistical Science 16:101-133 

Burgasser, A. J., Kirkpatrick, J. D., Reid, I. N., et al. 2003, ApJ, 586, 512 
Burgasser, A. J., Kirkpatrick, J. D., Cruz, K. L., et al. 2006, ApJS, 166, 585 
Burgasser, A. J., Reid, I. N., Siegler, N., et al. 2007, Protostars and Planets V, 427 
Burningham, B., Pinfield, D. J., Lucas, P. W., et al. 2010, MNRAS, 406, 1885 
Burningham, B., Cardoso, C. V., Smith, L., et al. 2013, MNRAS, 433, 457 
Casewell, S. L., Dobbie, P. D., Hodgkin, S. T., et al. 2007, MNRAS, 378, 1131 
Chiu, K., Fan, X., Leggett, S. K., et al. 2006, AJ, 131, 2722 
Cushing, M. C., Kirkpatrick, J. D., Gelino, C. R., et al. 2011, ApJ, 743, 50 
Close, L. M., Zuckerman, B., Song, I., et al. 2007, ApJ, 660, 1492 
Delfosse, X., Tinney, C. G., Forveille, T., et al. 1997, A&A, 327, L25 
Dobbie, P. D., Kenyon, F., Jameson, R. F., et al. 2002, MNRAS, 335, 687 
Duchene, G., Bouvier, J., Moraux, E., et al. 2013, A&A, 555, A137 
Dupuy, T. J., & Liu, M. C. 2011, ApJ, 733, 122 

Epchtein, N., de Batz, B., Capoani, L., et al. 1997, The Messenger, 87, 27 
Gelino, C. R., Kirkpatrick, J. D., Cushing, M. C., et al. 2011, AJ, 142, 57 


-18- 


Jameson, R. F., Dobbie, P. D., Hodgkin, S. T., &; Pinfield, D. J. 2002, MNRAS, 335, 853 

King, R. R., Parker, R. J., Patience, J., & Goodwin, S. P. 2012, MNRAS, 421, 2025 

Konopacky, Q. M., Ghez, A. M., Rice, E. L., & Duchene, G. 2007, ApJ, 663, 394 

Kraus, A. L., White, R. J., & Hillenbrand, L. A. 2005, ApJ, 633, 452 

Kraus, A. L., White, R. J., & Hillenbrand, L. A. 2006, ApJ, 649, 306 

Kraus, A. L., & Hillenbrand, L. A. 2007, AJ, 134, 2340 

Kraus, A. L., & Hillenbrand, L. A. 2012, ApJ, 757, 141 

Krist, J. E., Hook, R. N., & Stoehr, F. 2011, Proc. SPIE, 8127 

Lafreniere, D., Jayawardhana, R., Brandeker, A., Ahmic, M., & van Kerkwijk, M. H. 2008, ApJ, 
683, 844 

Luhman, K. L. 2004, ApJ, 602, 816 
Luhman, K. L. 2007, ApJS, 173, 104 

Luhman, K. L., Mamajek, E. E., Allen, P. R., Muench, A. A., & Finkbeiner, D. P. 2009, ApJ, 691, 
1265 

Luhman, K. L., Allen, P. R., Espaillat, C., Hartmann, L., & Calvet, N. 2010, ApJS, 186, 111 
Luhman, K. L. 2012, ARA&A, 50, 65 

Lodieu, N., Deacon, N. R., & Hambly, N. C. 2012, MNRAS, 422, 1495 

Liu, M. C., Leggett, S. K., Golimowski, D. A., et al. 2006, ApJ, 647, 1393 

Liu, M. C., Dupuy, T. J., & Leggett, S. K. 2010, ApJ, 722, 311 

Liu, M. C., Deacon, N. R., Magnier, E. A., et al. 2011, ApJ, 740, LL32 

Liu, M. C., Dupuy, T. J., Bowler, B. P., Leggett, S. K., & Best, W. M. J. 2012, ApJ, 758, 57 

Melis, C., Reid, M. J., Mioduszewski, A. J., Stauffer, J. R., & Bower, G. C. 2014, arXiv: 1408.6544 

Martin, E. L., Basri, G., Zapatero-Osorio, M. R., Rebolo, R., & Lopez, R. J. G. 1998, ApJ, 507, 
L41 

Martin, E. L., Brandner, W., Bouvier, J., et al. 2000, ApJ, 543, 299 

Martin, E. L., Barrado y Navascues, D., Baraffe, I., Bouy, H., & Dahrn, S. 2003, ApJ, 594, 525 
Moraux, E., Bouvier, J., Stauffer, J. R., & Cuillandre, J.-C. 2003, A&A, 400, 891 


-19- 


Nagashima, C., Dobbie, P. D., Nagayama, T., et al. 2003, MNRAS, 343, 1263 
Neuhauser, R., Brandner, W., Alves, J., Joergens, V., & Comeron, F. 2002, A&A, 384, 999 
Pecaut, M. J., Mamajek, E. E., & Bubar, E. J. 2012, ApJ, 746, 154 
Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9 

Press, W. H., Teukolsky, S. A., Vetterling, W. T. and Flannery, B. P., 1988, Numerical Recipes in 
C, Second Edition, Cambridge University Press, NY 

Reid, I. N., Kirkpatrick, J. D., Gizis, J. E., et al. 2000, AJ, 119, 369 

Reid, I. N., Cruz, K. L., Kirkpatrick, J. D., et al. 2008, AJ, 136, 1290 

Reid, I. N., & Hawley, S. L. 2000, New light on dark stars. Red dwarfs, low-mass stars, brown 
dwarfs., by Reid, I. N.; Hawley, S. L.. Springer, London (UK), 2000, ISBN 1-85233-100-3, 
Published in association with Praxis Publishing, Chichester (UK). 

Radigan, J., Jayawardhana, R., Lafreniere, D., et al. 2013, ApJ, 778, 36 

Stauffer, J. R., Schultz, G., & Kirkpatrick, J. D. 1998, ApJ, 499, L199 

Stauffer, J. R., Hartmann, L. W., Fazio, G. G., et al. 2007, ApJS, 172, 663 

Todorov, K., Luhman, K. L., & McLeod, K. K. 2010, ApJ, 714, L84 

Todorov, K. O., Luhman, K. L., Konopacky, Q. M., et al. 2014, ApJ, 788, 40 

van Dokkum, P. G. 2001, PASP, 113, 1420 

Zapatero Osorio, M. R., Galvez Ortiz, M. C., Bihain, G., et al. 2014, A&A, 568, AA77 
Zapatero Osorio, M. R., Bejar, V. J. S., Martin, E. L., et al. 2014, arXiv:1410.2383 


This preprint was prepared with the AAS IAT^X macros v5.2. 



- 20 - 


F814W 



F850LP 



TinyTim 












1 

1 














Semi-Empirical 


o.oio 

0.005 

0.000 


-0.005 


- 0.010 

-0.015 



Fig. 1.— The average residuals of all WFC3 images after fitting the original Tiny Tim model (top) 
and fitting our semi-empirical PSF model (bottom). For viewing purposes, we display the average 
residuals as the normalized “(image-model)” in each filter. When using our original Tiny Tim 
model, the average residuals results in concentrated residual flux at instrumental position angles of 
30 — 50° and 150 — 170° for both the F814W and F850LP filters. This would bias our binary fitting 
technique to preferentially recover companions with these position angles. Therefore we computed 
a semi-empirical model PSF using the original Tiny Tim model as a starting point as detailed in 
T3.ll The resulting average residuals are improved by a factor of 4-5x from ~ 5% and ~ 6% to 
~ 0.9% and ~ 2.3% in F814W and F850LP respectively. The residuals are also smoother, no longer 
containing concentrations at position angles of 30 — 50° and 150 — 170°. 
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Fig. 2.— Number density of binary parameters returned when fitting images of single stars with 
a binary PSF model, i.e., false positive detections, for the median S/N F814W case. Overall, the 
most common false positives have wide with separations of >3 pixels, and faint flux ratios of >4 
mag companions, but at separations of <1 pixels (<0 / /04), the majority of false positives range with 
flux ratios of 0 — 3 mag. The dotted red histograms are the initial guesses for the false positives 
uniformly distributed in log separation and flux ratio. The 1% false positive curve ( H3.4[) is over 
plotted (diamonds). 
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Fig. 3.— Illustration of our calculation of the false positive curve using the case of artificial 
binaries at separations 0.79 — 1.0 pixels (« 0 / . / 03 — 0 / . , 04) as an example. TOP LEFT: The white 
histogram is the distribution of the recovered flux ratios Am for artificial binaries at separations 
of 0.79 — 1.0 pixels. The histogram with slashes are false positives recovered by using our binary 
fitting technique on single star images. The vertical red dashed line is the median false positive 
flux ratio. MIDDLE LEFT: We normalize the histogram of false positive flux ratios (slashes) to 
the white histogram of recovered artificial binary flux ratios by conservatively assuming that any 
artificial binaries with recovered flux ratios larger than the median false positive flux ratio (vertical 
dashed red line) are most likely false positives themselves. BOTTOM LEFT: 1% (black solid 
line), 5% (red) and 10% (blue) false positive fractions as a function of flux ratio. RIGHT: The 
false positive curve is constructed by repeating the process for all separation bins. The stars denote 
the Am corresponding to 1%, 5% and 10% false positive fraction at separations of 0.79 — 1.0 pixels 
shown at the bottom left. 
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Fig. 4.— False positive curves computed at the minimum (triangles) and median (squares) signal- 
to-noise of our WFC3 images of Pleiades brown dwarfs ( 1 13.41) for the F814W (black) and F850LP 
(grey) filters. As expected the minimum S/N false positive curves have brighter false positives than 
the median S/N curves in a given filter. 
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Fig. 5.— Contrast curves for the median F814W S/N case constructed using our default and a 
stricter companion recovery criteria (see 1)3.51) . In our default criteria, we required companions be 
recovered to within 1 mag of the input binary parameters (dark grey solid line). We can recover 
companions with flux ratios <1 mag and separations >0"04. We also tested a stricter criteria, 
and required recovered companions to be within within 0.3 mag of the input flux ratio (black solid 
line). The contrast curves are identical for separations >0 / /055. With the stricter recovery criteria, 
companions with separations <0 / /04 and flux ratios < 0.5 mag were detectable. Both contrast 
curves required that recovered artificial binaries be within 0.2 pixels of input (x, y ) position. We 
adopt our default criteria given that most brown dwarf binaries are found to have near unity flux 
ratios. The 1% false positive curve is shown for comparison (light grey dotted line). The contrast 
curve drops at 200 mas due to difficulty in fitting artificial binaries at the edge of our cut-out region 
of ±10 pixels. 
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Fig. 6.— 90% completeness contrast curves for our F814W and F850LP observations of 11 young 
L dwarfs (<40 M Jup ) in the Pleiades. Our contrast curves rule out the majority of expected brown 
dwar f binaries, given that most binaries in the field have mass ratios >0.6 and separations <25 
AU (lBurgasser e t al.l 20071 1. We convert our detection limit fl ux ratios in WFC3 bandpasses to 


mass ratios using the dis tance to the Pleiades (13 6.2 pc. 


from lBaraffe et al.l ( 20031 ) tied to BT-Settl models ( Allardl 120141 ) 


Melis et all 12014] ) and evolution models 
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Fig. 7.— Example completeness maps for several our 11 young L dwarfs (<40 M Jup ) in the Pleiades. 
At each point, the completeness map represents the percentage of binaries that would have been 
recovered given our observations. We convert our detection limit flux ratio s in WFC3 bandpasses 
to m ass ratios using the d istance to the Pleiades 113 6.2 nc.lMelis et al.ll2014h . and evolution models 
from Baraffe et al. ( 2003 ) tied to BT-Settl models ( Allard 2014). 
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Fig. 8.— The wide (>10 AU) brown dwarf binary frequency as function of the age for young 
star forming regions, the intermediate age Pleiades, and the held (see ( 14.2)1 . All populations are 
shown for a common mass range (25 — 40 M Jup ). Low mass brown dwarf binaries may very well be 
infrequent across a wide range of ages. 












